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The effects of nonmagnetic impurities on the neutron scattering intensity are studied for a model of 
the copper oxide layers in the normal state. The contribution to the Q = (n, n) neutron scattering 
intensity from processes involving the scattering of the spin-fluctuations from an impurity with large 
momentum transfers are calculated within the random phase approximation. It is shown that this 
type of scatterings could lead to a peak in the neutron scattering intensity in the normal state. 

PACS Numbers: 74.62. Dh, 74.72.Bk, 75.10.Lp, 74.25.-q 



An important feature found by the neutron scattering 
experiments on YBa2Cu307 is a resonant peak which is 
observed at 41 meV energy and Q = {%, 7r) momentum 
transfers Jl] ||. In this material, the peak is observed 
only in the superconducting state. The theories which 
have been put forward attribute the resonant peak to an 
instability in the particle-particle [||J|] or the magnetic 
channels It is of interest to determine how much 

each of these two channels contribute to the resonant 
peak in YBa 2 Cu30y. 

Nonmagnetic Zn impurities have been used as a probe 
of the magnetic correlations in NMR experiments, and 
yielded valuable information Neutron scattering ex- 
periments on Zn substituted YBa2Cu30y also find that 
the Zn impurities have important effects on the spin fluc- 
tuation spectrum [jlO 11|. Upon substituting 0.5% Zn, 
the resonant peak of the superconducting state becomes 
observable in the normal state as a broadened peak cen- 



tered at 40 meV with a width of 10 meV 11 1. As the tem- 
perature is lowered, the peak intensity increases smoothly 
without any observable discontinuity in the slope at the 
superconducting transition at T c — 87°K. It is found that 
significant amount of spectral weight develops in the peak 
already in the normal state. For instance, at 94°K the 
amplitude of the neutron scattering intensity is about 
half of its value at 22°K. In this sample, there is no ob- 
servable spectral weight at frequencies less than ~ 35 
meV. However, when the concentration of the impuri- 
ties is increased to 2%, the neutron scattering intensity 
develops weight at frequencies as low as 5 meV, and a 
broad peak is observed at about 35 meV 10 1. The mo- 
mentum width of the structure at 35 meV is about 0.5 
A -1 , which is twice that of the peak in the sample with 
0.5% Zn. These are the main features of the changes in- 
duced by the Zn impurities. Understanding the reason 
for why dilute Zn impurities induce a peak in the neutron 
scattering intensity in the normal state would be useful 
for understanding the origin of the resonant peak in pure 
YBa 2 Cu 3 7 [0. 

Here, the effects of dilute Zn impurities on the neutron 
scattering spectral weight Imx(Q, u>) will be studied for 
a two-dimensional model of the layered cuprates in the 
normal state. The contributions to Imx(Q = (tt,it),u>) 



from processes involving the scattering of the spin fluctu- 
ations from an impurity with large momentum transfers 
near 2kp , where kp is the Fermi momentum, are cal- 
culated within the random phase approximation (RPA). 
It will be seen that this type of scatterings could lead 
to a peak in Im^(Q = (jt,tt),lli) in the normal state at 
loq = 2 1 /i |, where fj, is the chemical potential. 

During the scattering of the spin fluctuations from an 
impurity, momentum is not conserved and, hence, the 
Q = (n, 7r) component of the spin fluctuations could get 
mixed with the other wave vector components. For in- 
stance, if momentum Q* is transfered during a scattering, 
then the Q = (jr, n) component of the spin fluctuations 
would get mixed with the q = Q — Q* component. This 
process could be considered as an umklapp scattering of 
the spin fluctuations by the impurity potential. It has 
been previously noted that this type of scatterings with 
large momentum transfers could play a role in enhancing 
the uniform susceptibility as antiferromagnetic correla- 
tions develop in Zn substituted YBa2Cus07_5 [Q. It 
has been also shown that without the umklapp scatter- 
ings, the main effect of the impuritites on Im%(Q, uS) is 
to cause a smooth smearing Here, it will be shown 
that, when included within RPA, the scattering of the 
spin fluctuations with momentum transfers Q* ~ 2kp 
can lead to a peak in Imx(Q, ui) near loq in the normal 
state. The underlying reason for this is that there is a 
kinematic constraint against creating a particle-hole pair 
with center-of-mass momentum q = (tt, tt) — 2k^? and en- 
ergy u> > ujq. On the other hand, the scattering of the 
spin fluctuations with momentum transfers away from 
2kp has a smaller effect on Imx(Q,t^ ~ ^o)- 

In the following, the two-dimensional Hubbard model 
will be used to model the magnetic correlations of the 
pure system, and a static potential will be used to model 
the effective interaction between the electrons and a non- 
magnetic impurity. The two-dimensional single-band 
Hubbard model is defined by 
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where Cj CT (ct ) annihilates (creates) an electron with spin 
a at site i, i is the near-neighbor hopping matrix element 
and U is the onsite Coulomb repulsion. For simplicity, 
the hopping t, the lattice constant a and h are set to 1. 

When an impurity is introduced at site ro of the sys- 
tem described by Eq. (Q), the translational invariance is 
broken and, in this case, the magnetic susceptibility is 
defined as 

X (q,q', iWm ) = / dre^" T (m-(q,r)™ + (q',0)), (2) 
Jo 

where m+(q) = ^ _1/2 E P c p+ q T c pl' m ~( t l) = 
(m + (q))T, and uj m = 2mnT. The susceptibility defined 
by Eq. (|^) has both diagonal q = q' and off-diagonal 
q ^ q' terms. If a static effective interaction is assumed 
between the impurity at site ro and the electrons, then 
the RPA expression for x m th e Hubbard model becomes 

X(q,q',iw m ) =xo(q,q',^m) 

+U^2xo(q,<l",iu m )x(q.",<l',iu m ), (3) 

q" 

where XoCQiQ' >^m) is the irreducible susceptibility 
dressed with the impurity scatterings, but not with the 
Coulomb correlations. Here, the diagonal q = q' terms 
of the irreducible susceptibility Xo(q, q', iw m ) will be ap- 
proximated by the Lindhard function of the pure sys- 
tem, xd i^m). This is a good approximation since 
it is already known from previous calculations that the 
effect of impurity scattering on the diagonal irreducible 
susceptibility is to cause a smooth suppression ||l3| , and 
this has a small effect compared to the nearly singular 
contribution arising from the umklapp processes. The 
q/q' off-diagonal terms of Xo(q, q', iu m ) will be calcu- 
lated in the lowest order in the strength of the effective 
electron-impurity interaction as shown diagrammatically 
in Fig. 1(a). It will be seen that, within this model, 
the impurity induced changes in the Q = (it, tt) spin- 
fluctuation spectrum are closely related to the spectrum 
of the irreducible off-diagonal susceptibility xo(Q, q, w) 
with q = Q — 2kp. After including the effects of the 
impurity scatterings and the Coulomb correlations, the 
impurity averaging is done by averaging over the impu- 
rity location ro, which sets q = q'. This way, one obtains 
x(Q, Q, iu m ), which gives the neutron scattering spectral 
weight after analytic continuation iu) m — > u> + iS. In the 
following, x(Qi Q,iu m ) will be denoted by x(Q,«w m ). 

A static effective electron-impurity interaction for an 
impurity located at site ro can be written as 

k p,cr 

which involves scatterings with all possible wave vectors. 
It has been noted that the electronic correlations play a 
role in determining the structure of Vk Jl^| . It will be 
seen that for u « wq the important scatterings are the 
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FIG. 1. (a) Feynman diagrams representing the low- 
est order terms of the irreducible umklapp susceptibility 
Xo(Q,q, u). (b) Sketch of a single-particle scattering pro- 
cess involving Q* = 2k^ momentum transfer in the Brillouin 
zone. Here ki? is the Fermi wave vector along (1, 1) and the 
solid curve represents the Fermi surface of the noninteracting 
system at finite hole doping. 

ones which transfer Q* = 2kp, since in this case 
Xo(Q,q, w) with q = Q Q* has nearly singular be- 
haviour. For other values of Q*, Xo(Q> <F w ) is a smooth 
function of u> with small amplitude. Hence, Eq. (|3|) will 
be solved by using only the Q* component of Eq. ([|), 

v °J2 c l+Q*° c *><» ( 5 ) 

per 

where Vq = Vq» is taken as a parameter. Since the 
important umklapp contributions originate from a small 
range of values of Q* near 2kp, this is a reasonable ap- 
proximation. This procedure for solving Eq. (||) was nec- 
essary, because the calculation needs to be carried out 
on a large lattice in order to have sufficient frequency 
resolution. A single-particle scattering process given by 
Eq. (J|) is illustrated in Fig. 1(b) for Q* = 2kp momen- 
tum transfer in the Brillouin zone. 

Within the presence of the scattering term given by 
Eq. (||) , the solution of the RPA expression, Eq. (||) , be- 
comes 

x(Q,w) - d- 1 |xo : (Q,w)(i - tfxo(q.w)) 

+4C/( X o(Q ! q,^)) 2 |, (6) 

where D = (1 - Uxo(Q,w))(l - Uxfifau)) - 
4(t^xo(Q,q, w)) 2 and q = Q - Q*. Here, the factor 
of 4 multiplying (xo(Q, q, u)) 2 is to take into account 
the scatterings with momentum transfers (±Q*,±Q*) 
and (±Q*,=fQ*), where Q* = (Q*,Q*). The expression 
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for the irreducible umklapp susceptibility Xo(QjQ: ik>m), 
which is illustrated in Fig. 1(a), is given by 



Xo 



(Q,q_,iu m ) = -V — ^2 |G (p + Q,iw„ + ?w m ) 

p,iuj n 

xG (p, iu; n )G (p + Q q, iLu n ) 
+G (p, iu> n )G (p + q, iuJ n + tUm) 

xG a (p+Q,iLu n + iuj m )\. (7) 



Here, Go(p, iu> n ) — (iui n — e P ) _1 with u) n = (2n + 1)ttT 
and £ p = —2t(cosp x + cosp y ) — jjL is the single-particle 
Green's function of the pure system. After summing over 
iuj n in Eq. (j^) analytically, the resulting sum over p is 
carried out numerically. The Lindhard susceptibility Xo 
used in Eq. <M) is given by 



xf(k,w) 



1 



N ^ w - (e F 



/(£ P +k) - /(e P ) 



-k - £ P ) + 



(8) 



which is also evaluated numerically. The replacement of 
the diagonal irreducible susceptibility by the Lindhard 
susceptibility in solving Eq. (||) does not affect the re- 
sults presented here. It has been shown in Ref. |l3| that 
scatterings from an onsite impurity potential lead to a 
smooth smearing of the structure in the diagonal irre- 
ducible susceptibility. A similar effect was found for the 
u> = case in Ref. [fl2f using an extended impurity po- 
tential. 

The following results have been obtained by using 
U = 1.7, temperature T = 0.05 and filling (n) = 0.86. 
For these values of T and (n), one has luq = 2\n\ ks 0.55. 
Figure 2 shows results obtained for momentum transfer 
Q* = 2ki?, where kp was taken along (1,1) for simplic- 
ity. In Fig. 2(a), the real and the imaginary parts of 
Xo(Q,q, w) are shown. For comparison, in Fig. 2(b), the 
u> dependence of the Lindhard susceptibility Xo (Q,<^) is 
shown. The solid line in Fig. 2(c) represents Imx(Q,w) 
for Vo — 0.05. Also shown in this figure by the dashed 
line is the result for Vo = 0, which corresponds to the 
pure case. Here, it is seen that Imxo(Qiq 5 w ) is finite 
for lu < loq with a sharp cut-off at luq. The cut-off in 
Imxo(Qj q> w ) is due to the fact that there is a kinematic 
constraint against creating a particle-hole pair with the 
center-of-mass momentum q = {%, n) — 2kp and energy 
lu > loq. In the corresponding real part there is a sharp 
peak at loq. Both the real and the imaginary parts of 
Xo(Q,q, w) play a role in suppressing the denominator 
of the expression in Eq. (](]), and in leading to the peak 
in Imx(Q,w) at loq- In this model, the peak is origi- 
nating from the magnetic channel, since it is due to an 
instability of Eq. (||) . 

The spectrum of Xo(Q,q, w ) seen in Fig. 2(a) is 
quite different than that of the Lindhard susceptibil- 
ity xo(Q, w ) seen in Fig. 2(b). In the limit T — > 0, 
Imxo(Q: w ) develops a gap for to < loq. On the other 
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FIG. 2. (a) Irreducible umklapp susceptibility Xo(Q, q, u>) 
versus u normalized by Vo. Here q = Q — Q*. (b) Lindhard 
susceptibility Xo(Qi w ) versus to. (c) Spin-fluctuation spectral 
weight Imx(Q,^) versus ui for Vo = 0.05 and 0.0. These 
results were obtained for Q* = 2k^. 

hand, Imx (Qi <3i w ) is finite only for lu < loq. Hence, the 
umklapp scatterings could contribute to Imx(Q,w) for 
lu < u)q. The quantitative features of Imx(Q,w) depend 
on the model parameters. For instance, if U is increased 
from 1.7 to 2.0, the hump in Imx(Q,cj) centered at lu — 
0.5 grows and the dip below the peak gets smaller. The 
scattering of the quasiparticlcs by the spin fluctuations 
would also affect the quantitative features. 

Figures 3(a) and (b) show the temperature evolution 
OI Xo(Q ; q7 w ) an d Imx(Q,w) for V = 0.05 on a smaller 
lu scale. As T increases, the peak in Imx(Q,w) disap- 
pears. Figure 3(c) shows the development of the peak in 
Imx(Q,w) as Vo is turned on at T = 0.05. 

In order to have further information on the momen- 
tum and frequency structure of xo(Q,q7 ti; )j results on 
its real part are shown in Fig. 4. Figure 4(a) shows 
Rexo(Q,q, oj) versus lu obtained at values of Q* which 
are slightly different than 2kp, but still along (1,1). Here, 
the numbers next to the curves indicate the magnitude 
of Q* in units of 2k p. As Q* is varied away from 
2ki?, the resonance frequency shifts and the amplitude 
of Rexo(Q,q7 w ) decreases. Hence, when all the contri- 
butions from scatterings with momentum transfers near 
2ki? are taken into account, one expects the peak in 
Imx(Q;k>) at luq to have a finite width. However, it is 
not possible to estimate the total spectral weight in the 
peak. For instance, the detailed k dependence of Vk for 
k ~ 2k^ is not known. 
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FIG. 3. Temperature variation of (a) Xo(Q> 1, w)/Vo and 
(b) Imx(Q,^) versus lo for Vo = 0.05. In these two figures, 
the chemical potential has been kept fixed, while increasing 
the temperature, (c) ImxCQjW) versus lo for various values 
of Vo at T — 0.05. These results were obtained for Q* = 21cf. 
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FIG. 4. Re xo(Q, q, w)/Vo versus lo (a) for Q* near 2kF 
and (b) for various values of Q* away from 2kj? at (n) = 0.86. 
(c) Rexo(Q, q,u>)/Vo versus lo for Q* near 2ki? at (n) = 0.7. 



In Fig. 4(b), Rexo(Q,q,^) is shown for various val- 
ues of Q* away from 2kp. Here one observes that 
R- e Xo(Qi w ) is a smooth function of lu and has a much 
smaller amplitude near loq than in the case for Q* rs 2hp. 
Other values for Q* away from 2kp were also tried, but 
results similar to those in Fig. 4(b) were obtained. This 
supports keeping only the Q* « 2kp component, Eq. (||), 
of the effective impurity interaction, Eq. (|J), in solving 
the RPA equation for xiQ, 1 ^) at w near cj - Finally, 
Fig. 4(c) shows Rexo(Q,q, w) versus cj for Q* near 2kp 
at filling (n) = 0.7, where the corresponding uj is 1.44. 
Comparing with the results shown in Fig. 4(a), one ob- 
serves that the qualitative features do not change with 
the band filling. 

The neutron scattering experiments by Fong et al. |Tl| l 
have found that, as the temperature is lowered below 
250°K, a peak develops at w = 40 meV in Imx(Q,^) 
of 0.5% Zn substituted YBa2Cu307. Here, the effects of 
dilute nonmagnetic impurities on Imx(Q,w) have been 
studied within a two-dimensional model in the normal 
state. It has been found that the impurity induced 
changes in Imx(Q, are closely related to the spectrum 
of the irreducible off-diagonal susceptibility xoiQ,^.,^)- 
The possibility of the structure in Xo(Q> CU k>) at w ~ loq 
leading to a peak in Imx(Q,w) has been pointed out. 
Within this model, the peak is due to the impurity scat- 
tering of the spin fluctuations with momentum trans- 
fers near 2kp. However, in obtaining these results the 
Coulomb correlations were treated within the random 
phase approximation, and only the lowest order terms 
of Xo(Qi Qj w ) have been included. It remains to be seen 
by how much the higher order corrections will modify the 
results discussed here. 

The author thanks P. Bourges, H.F. Fong and B. 
Keimer for helpful discussions. The numerical computa- 
tions reported in this paper were performed at the Center 
for Information Technology at Koc University. 
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